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SECTION  I  -  INTRODUCTION 


The  AFCWC  automated  cloud  analysis  model  has  been  operational  since 
January  1970.  This  model,  known  as  the  Three-Dimensional  Nephanalysis 
(3DNEPH),  relies  heavily  on  meteorological  satellites  to  provide  cloud 
information  over  the  entire  world.  Using  the  albedo  and  temperature  data 
inferred  from  reflected  light  and  emitted  infrared  imagery,  respectively,  the 
model  is  able  to  measure  accurately  the  amount  of  cloudiness  in  individual 
cloud  layers.  The  accuracy  of  this  technique  is  far  less  reliable  in  polar 
and  temperate  cone  winter  regions  than  elsewhere.  In  these  regions,  snow  and 
ice  fields  are  virtually  indistinguishable  from  low  level  clouds  in  both 
spectral  bands.  However,  several  investigators  suggested  that  near-infrared 
and  visible  band  sensors  might  be  used  to  improve  snow/cloud  discrimination. 

As  a  result,  the  Defense  Meteorological  Satellite  Program  (DMSP)  Program 
Office  funded  the  design  and  development  of  a  unique,  inexpensive  orbital 
instrument  to  exploit  this  new  technology.  We  were  designated  as  the  office 
of  primary  responsibility  for  developing  and  testing  new  snow/cloud 
discrimination  techniques  as  part  of  the  on-orbit  operation  of  the 
instrument.  Our  objectives  were  as  follows: 

Collect  aud  evaluate  data  from  the  DMPP  Snow/Cloud  Special  Sensor  (SSC)  to 
test  the  hypothesis  that  snow/cloud  scene  discrimination  is  possible  through 
interpretation  of  a  combination  of  SSC  and  DMSP  Operational  Linescan  System 
(0LS1  data. 

list  potential  benefits  if  cloud  and  snow  information  is  extracted  from 
S SC-type  imagery  and  processed  by  existing  models. 

We  discuss  the  characteristics  of  the  instruments  in  Section  2  and  provide 
a  detailed  description  of  our  data  evaluation  in  Section  3.  Section  4  reports 
the  factors  which  limited  the  utility  of  the  sensors.  We  list  our  conclusions 
and  the  potential  benefits  of  processing  near- infrared  satellite  data  in 
Sections  5  and  6,  respectively.  Section  7  is  a  list  of  references. 
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SECTION  2  -  INSTRUMENTATION 

2.1  OPERATIONAL  LINESCAN  SYSTEM  (OLS) 

The  OLS  is  a  duel  channel  scanning  radiometer  which  collects 
reflecced  light  (L)  ana  emitted  thermal  (T)  infrared  energy  in  the  0.4  to  l.i 
um  and  10.2to  12.8  um  spectral  bands,  respectively.  Energy  is  captured  by  a 
scanning  telescope  system  driven  in  a  sinusoidal  motion  by  counter-reacting 
coiled  springs  and  pulsed  motor.  The  incoming  beam  of  energy  is  split  into 
two  separate  paths  and  then  focused  onto  discrete  L-  and  T-detector  elements. 
The  system  stoves  the  instantaneous  field  of  view  (1FOV)  of  these  elements 
across  the  subtrack  with  maximum  scanning  velocity  at  nadir  and  reversals  at 
the  edge  of  scan.  The  angular  1F0V  is  changed  dynamically  at  selected 
portions  of  scan  in  order  to  maintain  near  constant  resolution. 

Electronic  gair.  state  settings  are  used  to  adjust  L-detector 
sensitivity  to  approximate  changes  in  solar  irraaiance  caused  by  variation  of 
solar  zenith  angle  within  the  swath.  The  L-detector  is  designed  to  produce  an 
output  signal  strength  directly  proportional  to  incoming  scene  radiance.  This 
output  is  digitized  into  one  of  128  possible  greyshade  values  (i.e.,  /  bits). 
These  greyshades  are  truncated  upon  receipt  at  AFGVC  to  one  of  64  possible 
greyshade  values  (6  bits). 

The  T-detector  is  designed  to  produce  an  output  signal  strength 
directly  proportional  to  the  equivalent  Black  Body  (BB)  temperature  contained 
within  the  IFOV.  However,  atmospheric  attenuation  reduces  incoming  radiance 
so  that  equivalent  BB  temperatures  are  less  than  temperatures  of  radiating 
surfaces.  Figure  2.1  contains  a  plot  of  average  ambient  air  temperature 
versus  matching  T-detector  output  for  clear  scenes.  The  detector  response  io 
digitized  into  one  of  128  possible  greyshade  values  (7  bits).  These 
greyshades  are  truncated  upon  receipt  at  AFGMC  to  one  of  64  possible  greyshade 
values  (6  bits). 

The  L-  and  T-detector  data  are  received  in  near  constant  1.5  am 
resolution  at  AFGWC  where  they  are  smoothed  spatially  to  a  near  constant  3.0 
run  resolution.  The  resultant  light  and  thermal  IK  band  data  are  copied  to 
7-track,  800  BPI  magnetic  tapes.  Beginning  October  1979,  L-  and  T-detector 
data  were  saved  from  one  orbital  revolution  each  day  for  this  study. 

2.2  SNOW/CLOUD  SPECIAL  SENSOR  (SSC) 

The  SSC  is  a  narrow  spectral  band,  staring  mode,  "push  broom" 
scanning  radiometer.  This  radiometer  is  a  48  element  linear  array  of 
germanium  photovoltaic  detector  elements  which  lies  in  the  image  plane  of  a 
40.2°  wide  angle  lens  designed  to  provide  constant  angular  field  cf  view 
(FOV)  per  array  element  (14.6  mrad  along  track,  9.61  mrad  cross  track  on  14. f 
mrad  centers).  The  instrument's  FOV  extends  from  subtrack  toward  the  sun  and 
covers  about  4C0  nm  on  the  earth's  surface.  Figure  2.2  shews  the  nominal 
detector  footprint. 
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T-Detectcr  Greychadu  — ^ 


Figure  2.1.  Plots  of  average  apparent  Black  Body  temperature  (T  )  from  clear 
radiating  surfaces  and  corresponding  average  ambient  air  temperatures  (Tc)  versus 
matching  T-detector  greyshade  values. 
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The  instrument  collects  reflected  solar  energy  in  the  1.51  to  1.63  urn 
spectral  band  for  input  signal  strength.  Sixteen  gain  spates,  spaced  at  2.4 
dB,  are  available  for  use  at  specified  progressed  times  along  the  orbital 
track.  They  are  used  to  compensate  approximately  for  variations  in  scene 
illumination  from  near  subsolar  to  terminator  regions.  Channels  are  sampled 
in  the  "push  broom”  mode  every  two  seconds  and  are  individually  digitized  to 
six  bits  (64  possible  greyshades).  These  data  are  recorded  on  the  spacecraft, 
relayed  to  DHSP  Control  and  Readout  Stations,  shipped  to  AFGWC,  reformatted 
awl  recorded  on  tape  with  standard  time  and  epheoeris  information.  Table  2.1 
provides  additional  sensor  design  information. 
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ITEM 


VALUE 


UNITS 


Number  of  Channels 

Spectral  Band 

Output  Data  Quantisation 

Data  Frame  Time 

Each  Channel  Time  Constant 

Multiplex  Time  per  Channel 

IFOV  each  Channel 

Power  Consumption 

Weight 


48 

- 

1.51  to  1.63 

pm 

6 

Bits 

2 

sec 

0.68 

sec 

10 

msec 

14.6  x  9.1 

mrsd 

10.86 

watts 

8.58 

lbs 

Table  2.1.  SSC  System  Parameters 


SECTION  3  -  DATA  EVALUATION 


W«  began  early  orbit  check-out  of  SSC  shortly  after  launch  of  DKSP  WX5539, 
flight  four  (F4),  on  6  June  1979*  Although  the  aenaor  vaa  in  good  working 
order,  Veatinghouae  Electric  Corporation  (WEC)  engineera  reported  that  the 
ambient  aenaor  temperature  waa  alightly  higher  than  deirad.  Thia  high 
temperature  had  negligible  impact  on  aenaor  perfonaance.  Based  on  our 
findings  and  a  thorough  evaluation  by  system  engineera.  Space  Diviaion  (SD) 
declared  the  SSC  operational  on  10  June  79. 

Shortly  thereafter,  ayatea  engineera  uploaded  a  pre-programmed  gain  atate 
profile  to  coupeneate  for  variations  in  acene  illumination  throughout  the 
orbital  path.  Aa  a  reault,  the  detectora  produced  full  acale  output 
(aaturation)  whenever  acene  reflectances  (assuming  Lambertian  surfaces) 
exceeded  23Z. 

When  we  coapared  printer  plots  of  SSC  data  to  transparencies  of  matching 
1.5  na  resolution  DMSP  light  data,  we  noticed  that  the  sensor  usually 
saturated  if  liquid  water  clouds  filled  the  FOV.  The  SSC  output  froa 
snow-filled  FOVs  indicated  that  reflectivities  were  6X.  The  dynamic  range  of 
the  sensor  was  sufficiently  wide  to  handle  cloud-free,  snow-free  scene 
radiances  over  a  wide  variety  of  geographic  features.  However,  the  detector 
was  saturating  over  both  arid  and  semi-arid  scenes.  Ve  felt  that  liquid  water 
cloud  reflectivities  and  arid  scene  reflectivities  could  be  equivalent  under 
certain  conditions.  However,  we  could  not  study  this  aspect  without  expanding 
the  dynamic  range  of  the  aensor.  If  these  reflectivities  were  nearly 
equivalent  then  the  task  of  developing  reliable  cloud/no  cloud  decision 
algorithms  for  worldwide  use  could  be  extremely  difficult.  We  decided  to 
identify  areas  in  the  world  where  reflectivities  were  a  potential  problem  and 
to  develop  a  new  gain  state  profile  to  expand  the  dynamic  range. 

To  observe  worldwide  cloud-free  reflectivities,  we  considered  developing 
software  to  determine  where  cloud-free  SSC  data  were  acquired.  Ve  also 
considered  using  nephsnalyses  produced  by  the  5DNEPH  model  to  make  these 
determinations.  Both  of  these  approaches  were  ruled  out  because  the  manpower 
and  computer  costs  would  have  been  extremely  high.  We  could  get  comparable 
results  at  a  much  smaller  cost  by  computing  minimum  SSC  brightness  worldwide. 
The  SSC  minimum  brightness  is  equal  to  detector  response  from  cloud-free 
features,  except  in  those  cases  where  cirrus  clouds  are  less  reflective  than 
arid  terrain  features.  Ve  developed  software  to  compute  minimum  SSC 
brightness  on  an  AFGVC  eighth-mesh  scale  (25  nm  grid  spacing  at  60°N).  Ve 
used  July  and  August  1979  SSC  data  to  prepare  a  SSC  minimum  brightness  data 
base  and  displayed  the  results  on  a  photographic  transparency  (Figure  3.1). 

SSC  digital  responses  were  represented  by  16  possible  greyshade  values,  where 
black  and  white  correspond  to  minimum  and  maximum  responses,  respectively.  We 
identified  areas  in  North  Africa,  the  Middle  East,  Soviet  Union,  Mexico,  and 
the  United  States  where  reflectivities  were  sufficiently  high  to  cause  the 
sensor  to  saturate. 


Figure  3.1.  A  display  of  SSC  minimum  brightness  (next  page)  based  on  SSC  data 
collected  during  July  and  August  1979.  Brightness  is  represented  by  16  grey- 
shade  values  ranging  from  black  to  vhits.  The  sensor  was  set  to  produce  full- 
scale  output  whenever  scene  reflectances  exceeded  approximately  23Z.  The  white 
disk  at  the  pole  is  merely  a  data  void  area.  Notice  that  scai-arid  and  arid 
terrain  features  produced  full-scale  output.  Snow-covered  Greenland  is  much 
darker  in  comparison  to  the  semi-arid  regions  by  virtue  of  the  low  refectivity 
of  snow  in  the  near- infrared  portion  of  the  spectrum. 
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During  September  and  October  1979,  we  developed  formulae  and  software  Co 
produce  a  new  gain  aCate  profile  Co  increaae  uniformly  Che  dynamic  range  of 
the  sensor  throughout  the  orbital  track.  The  dynamic  range  was  increased  by 
7.2  dB  in  late  October  1979,  raising  Che  threshold  for  sensor  saturation  from 
232  to  532  reflectivity.  Further  changes  were  not  possible  for  areas  within 
AO. 7  degrees  of  the  sub-solar  point  by  virtue  of  sensor  design  limitations. 
Unfortunately,  most  of  the  highly  reflective  features  fell  within  those 
areas.  We  soon  found  that  the  increased  dynamic  range  did  not  completely 
eliminate  the  saturation  problem. 

To  isolate  areas  of  persistant  high  reflectivity,  we  prepared  a  SSC 
minimum  brightness  data  base  using  data  acquired  on  the  following  days: 

31  Oct  1979; 

1,  2,  3,  13,  14,  15,  16,  20,  21  and  30  Nov  1979; 

1,  2.  5,  6,  7,  8,  9,  12,  13,  14,  15,  16  and  17  Dec  1979. 

We  displayed  this  data  base  in  transparency  format  (Figure  3.2)  and  noticed 
that  arid  areas,  such  as  the  Sahara,  had  reflectivities  greater  than  532.  We 
also  noticed  when  comparing  raw  L-detactor  data  with  SSC  data,  that  some 
liquid  water  clouds  also  had  apparent  reflectivities  greater  chan  S3Z.  At 
that  point  we  ruled  out  worldwide  use  of  the  SSC  detector  as  a  "stand-alone" 
cloud  detection  device.  We  concentrated  on  developing  a  multi-spectral 
algorithm  which  would  use  SSC  data  in  conjunction  with  DMSP  OLS  data  to  fluke 
cloud/no  cloud  and  snow/no  snow  decisions. 

We  prepared  statistical  analyses  of  SSC  minimum  brightness  within  several 
areas,  land  only,  scattered  throughout  the  Northern  Hemisphere  (Figure  3.3). 

We  eliminated  every  point  north  of  69°N  from  consideration  because,  by  the 
end  of  the  period,  all  of  those  points  were  either  in  twilight  or  in  total 
darkness.  We  found  that  areas  containing  persistent  snow  cover  had 
reflectivities,  in  the  near-lR,  between  6  and  132.  We  also  matched  the 
minimum  brightness  display  (Figure  3.2)  with  weekly  Snow  and  Ice  Boundary 
Charts  (Figures  3.4  thru  3.8)  prepared  by  NCAA/ NESS.  Notice  the  excellent 
agreement  between  an  area  of  bright  snow  in  the  mid-western  part  of  the  United 
States  (Figure  3.6)  and  a  similar  area  of  low  near-IR  reflectivity  in  the  SSC 
minimum  brightness  display  (Figure  3.2).  These  comparisons  suggest  that  low 
|  neer-IR  reflectivites  alone  are  an  excellent  indicator  of  snow  cover.  This 
conclusion,  however,  is  not  valid  over  water  where  water  and  ice 
|  reflectivities  are  comparable. 

1 

S  Our  next  task  centered  on  developing  an  understanding  of  the  functional 
\  relationship  between  SSC  and  primary  detectors  (both  L  and  T)  responses  over 
\  cloud,  cloud-free,  and  snow  covered  scenes.  We  developed  software  to  compute 
the  average  L-  and  T-detector  responses  within  the  IFOV  of  each  SSC  detector. 
IWe  analysed  21  quarter  orbits  (16  N.  Hemisphere;  5  S.  Hemisphere)  of  SSC, 
jL-detector,  and  T-detector  data  and  concluded: 

Thick  cirrus  clouds  reflect  from  132  to  262  of  incident  solar  energy  in 
the  SSC  spectral  band. 
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Figure  3.2.  A  display  of  SSC  minimum  brightness  (next  page)  based  on  SSC  data 
collected  during  October,  November,  and  December  1979.  Brightness  is  repre¬ 
sented  by  16  greyshada  values  ranging  from  black  to  white.  The  sensor  was  set 
to  produce  full-  scale  output  whenever  scene  reflectances  acceded  approximate¬ 
ly  53X.  The  white  disk  at  the  pole  is  merely  a  data  void  area.  Notice  that 
aeai-arld  and  arid  features  caused  full-scale  output.  However,  the  spatial 
extent  of  this  output  is  much  leas  than  before  (Figure  3.1).  The  dark  areas 
over  land  have  reflectivities  roughly  equal  to  the  reflectivity  of  snow  in  the 
near-infrared  portion  of  the  spectrum  (Valovcin  (1978)).  In  fact,  these  areas 
represent  the  spatial  extent  of  the  Northern  Hemisphere  snow  cover  exception¬ 
ally  well  (Figures  3.4  through  3.8). 
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Figure  3.3.  Statistical  analyses  of  SSC  minimum  brightness  were  prepared  within 
the  areas  containing  circled  numbers. 
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Clear  ocean  scenes  reflect  from  1  to  21  of  incident  solar  energy  in  the 
SSC  spectral  band. 

Liquid  water  clouds  reflect,  in  most  cases,  more  than  53X  of  incident 
solar  energy  in  the  SSC  spectral  band. 

The  SSC  response  is  directly  proportional  to  L-detector  response  over 
clear,  snow-free  features. 

Ve  still  needed  to  develop  a  functional  relationship  between  SSC  and 
L-detector  responses  from  snow  (ice)  and  cloud-filled  scenes.  We  planned  to 
use  3DNEPH  and  SNODEP  node  Is  for  cloud  and  snow  depth  information, 
respectively,  to  separate  matching  SSC,  L-detector,  and  T-detector  samples 
into  cloud-free,  snow-filled,  and  cloud-covered  categories.  We  then  planned 
to  prepare  scatter  diagrams  of  SSC  versus  both  L-  and  T-detector  responses  for 
each  of  these  categories  and  to  develop  the  required  multi-spectral  processing 
algorithm.  Furthermore,  we  planned  to  test  the  performance  of  the  algorithm 
against  independent  sets  of  merged  L-detector,  T-detector,  and  the  SSC  data 
and  compute  the  statistical  reliability  of  the  multi-spectral  processing 
technique.  However,  the  sudden  deactivation  of  the  SSC  on  29  December  19^9 
and  subsequent  cessation  of  all  planned  SSC  testing  forced  ua  to  revise  < ur 
plans.  This  revision  was  begun  by  reviewing  the  quality  of  the  3DNEPH  a  d 
SNODEP  analyses  which  were  on  hand.  Until  that  time,  we  had  been  qualit 
controlling  and  correcting  deficiencies  in  3DNEPH  analyses  of  F4  data  on iy  to 
meet  operational  requirements.  Also,  oust  of  the  saved  L-  and  T-detect<  r  data 
covered  the  eastern  portion  of  North  America  where  snow  cover  was  generally 
restricted  to  remote  areas  of  Canada  and  conventional  observations  of  s  iow 
depth  were  very  sparse.  Consequently,  the  SNODEP  model  relied  almost 
exclusively  on  climatology  as  a  source  of  snow  depth  information.  As  a  result 
the  SNODEP  analyses  in  those  areas  were  not  as  reliable  as  we  would  have  liked 
them  to  be.  Therefore,  ve  decided  to  use  manually  prepared  verification  data 
bases  in  lieu  of  the  3DNEPH  and  SNODEP  analyses. 

We  manually  produced  a  verification  data  base  using  transparencies  of  SSC 
and  matching  L-  and  T-detector  data.  We  selected  the  following  four  quarter 
orbits  of  L-detector,  T-detector,  and  SSC  data  for  verification: 

Data  revolution  2552  on  3  Dec  1979  over  North  America  (Figure  3.9). 

Data  revolution  2680  on  12  Dec  1979  over  North  America  (Figure  3.10). 

Data  revolution  2680  on  12  Dec  1979  over  Antarctica  (Figure  3.11). 

Data  revolution  2860  on  25  Dec  1979  over  Eurasia  (Figure  3.12). 

We  took  the  following  steps  to  prepare  the  verification  data: 

We  prepared  transparencies  of  the  L-  and  T-  detector  data  which 
included  timing  marks  every  2  seconds. 

We  adjusted  SSC  data  to  compensate  for  solar  zenith  angle  and  actual 
gain  state  thereby  normalizing  sensor  response  to  0°  solar  zenith  angle. 


Figure  3.9.  A  display  of  matching  DM3?  visual  (3.0nm  resolution)  and  SSC  lmag 
ery  (next  page),  left  and  right  respectively.  These  data  were  collected  over 
North  America  on  3  December  1979  during  orbital  revolution  2552.  In  the  SSC 
display,  clouds  are  much  brighter  than  clear,  snow-covered  terrain.  However, 
liquid  water  clouds  are  brighter  than  cirrus  clouds  in  the  SSC  Imagery.  Cloud 
free  water  scenes  appear  dark  In  both  displays.  The  black  vertical  lines  In 
the  SSC  Imagery  represent  gaps  in  the  SSC  data  in  the  cross-track  direction 
(Figure  2.2).  The  Image  processing  technique  has  removed  some  of  these  gaps. 
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Figure  3.10.  A  display  of  matching  K13P  visual  (3.0nsa  resolution)  and  SSC  data 
(next  page),  left  and  right  respectively.  These  data  were  collected  over  North 
America  on  12  December  1979  during  orbital  revolution  2680.  Areas  of  clear, 
snov-covered  terrain  can  ba  found  in  the  center-left  portion  of  the  visual  imag¬ 
ery.  The  snow  has  roughly  the  same  brightness  as  nearby  cloud  cover.  The  snow 
is  much  darker,  however,  than  the  clouds  in  the  SSC  display,  making  snow/cloud 
discrimination  easier  with  the  SSC  imagery  than  with  visual  data.  The  upper 
portion  of  the  SSC  imagery  (i.e.,  the  white  area)  did  not  receive  sufficient 
daylight  to  permit  image  processing.  All  but  the  upper  portion  of  the  display 
was  used  in  preparing  a  verification  data  base. 


Figure  3.11.  A  display  of  matching  DHSP  visual  (3.0na  resolution)  and  SSC  data 
(next  page),  left  and  right  respectively.  Theae  data  were  collected  over  Ant¬ 
arctica  on  12  December  1979  during  orbital  revolution  2680.  Notice  the  excel¬ 
lent  contrast  between  bright  clouds  and  dark  snow  in  the  SSC  display  compared 
to  the  relatively  poor  contrast  between  bright  snow  and  clouds  in  the  visual 
Imagery. 
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Figure  3.12.  A  display  of  matching  CMSP  visual  (3.0nm  resolution)  and  SSC  data 
(next  page),  left  and  right  respectively,  These  data  were  collected  over  Eur¬ 
asia  on  25  December  1979  during  orbital  revolution  2350.  A  narrow  band  of 
status  clouds  la  located  In  the  upper  quarter  of  the  visual  imagery.  Next  to 
these  clouds  Is  an  area  of  clear,  enow-covered  terrain.  BMSP  thermal  data  (not 
shown)  Indicate  that  the  cloud  top  temperature  la  warmer  than  the  adjoining  snow- 
covered  terrain.  Current  3DNEPH  satellite  data  processing  techniques  would 
have  failed  to  detect  these  clouds.  Notice,  however,  how  well  these  clouds 
stand  out  against  the  dark  snow  background  in  the  SSC  display.  Cirrus  clouds 
are  less  bright  than  liquid  water  clouds  in  the  SSC  imagery.  The  lover  portion 
of  the  images  cover  part  of  the  Middle  East.  Notice  that  clear,  arid  terrain 
features  caused  full-scale  output  of  -:he  SSC  sensor. 
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We  mapped  SSC  data  into  an  L-  and  T-detcctor  format  and  displayed  the 
results  on  transparencies. 

We  prepared  a  template  of  the  IFOV  of  the  SSC  detector  display  at 
nominal  altitude  to  match  L-  and  T-detector  displays. 

We  placed  the  template  over  the  L-  and  T-detector  'isplays  at  time 
marks  corresponding  to  known  SSC  data  times. 

We  inspected  all  available  data  and  decided  whether  there  was  snow, 
cloud,  or  no  cloud  within  the  IFOV  of  each  SSC  display.  We  treated  any  amount 
of  cloud  within  the  IFOV,  even  a  single  L-  or  T-detector  pixel,  es  a  cloudy 
SSC  pixel.  Similarly,  any  amount  of  anow  within  a  clear  IFOV  was  deemed  s 
snow-filled  pixel.  In  other  words,  we  made  no  provision  for  partly  cloudy 
scenes,  or,  for  that  matter,  mixed  snow  and  non-snow  scenes.  Mixed  snow  sad 
cloud  scenes  were  treated  aa  cloudy  scenes. 

We  used  these  results  to  stratify  the  merged  data  aeta  into  clear, 
cloudy,  and  snow  (ice)  categories. 

We  then  prepared  scatter  diagrams  of  responses  from  coincident  SSC  and 
average  L-detector  samples  and  coincident  SSC  and  average  T-detector  samples 
within  each  of  the  categories  (Figures  3.13  through  3.18). 

Analysis  of  these  diagrams  showed  the  best  potential  for  anow,  cloud,  and 
cloud-free  scene  discrimination  could  be  achieved  by  interpreting  coincident 
average  L-detector  and  SSC  samples.  Figure  3.13  reinforced  our  previous 
conclusion  that  cloud-free  SSC  response  is  directly  proportional  to  coincident 
average  L-detector  response.  Comparison  of  Figure  3.13  with  Figure  3.15 
showed  excellent  separation  between  clusters  of  points  of  SSC  versus  average 
L-detector  samples  acquired  over  cloud-free  and  snow-covered  scenes.  Detailed 
analysis  of  the  data  in  Figures  3.13,  3.15,  and  3.17  showed  that  clusters  of 
points  of  SSC  versus  average  L-detector  samples  acquired  over  clouds  were 
generally  located  away  from  the  clusters  of  clear  ar.d  snow-filled  points. 
However,  a  limited  amount  of  cloud  samples  overlap  slightly  into  the  scatter 
of  both  snow  and  cloud-free  samplea.  We  formed  a  composite  diagram  of  all 
samples  and  divided  the  matrix  of  possible  combinations  of  SSC  versus  average 
L-detector  responses  into  cloud,  snow  (ice),  and  cloud-free  regions  to 
minimize  the  overlap  (Figure  3.19).  The  matrix  provided  a  mechanism  for 
making  unique  cloud,  snow,  or  cloud-free  decisions  baaed  merely  on  matching 
SSC  and  average  L-detector  responses.  We  developed  an  automated  algorithm 
which  effectively  uses  the  matrix  a.«  a  sample  look-up  table  to  determine  scene 
content. 

During  the  final  phase  of  our  evaluation,  wa  determined  the  ctatistical 
reliability  of  this  algorithm  by  comparing  results  obtained  from  the  algorithm 
to  manually  prepared  verification  results.  Tables  3.1  through  3.4  contain 
results  obtained  from  the  algorithm  for  cloud-free,  snow-filled,  or 
cloud-filled  categories  versus  verified  results  for  each  of  the  four  selected 
quarter-orbits  of  data.  Table  3.5  depicts  the  combined  results.  Each  row 
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Figure  3.13.  Scatter  diagram  cf  L-detector  and  SSC  output  from  clear  scenes. 
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SCATTER  DIAGRAM  Of  RAW  F  D£  Q  OCCURRENCE  OF  T-DLTECTOR  GREY SHADES  fY-AXlSI 
VERSUS  SSC  GREYSHADES  (X-AXIS)  COLLECTED  OVER  CLEAR  SCENES 
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figure  3.15.  Scatter  diagram  of  L-detector  and  SSC  output  from  snow  scenes. 
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Figure  3.17.  Scatter  diagram  of  L-dotector  and  SSC  output  from  cloud  scenes. 


-33- 


r CATTER  DIAGRAM  OF  RAW  FREO  OCCURRENCE  OF  T-DETECTOR  GRL ^SHADES  (Y-AXIS) 
VERSUS  SSC  GPEYSHADES  (X-AXIS)  COLLECT FD  OVER  CLOUD  SCENES 


Q 

8 


I* 

1 

2 

7 

4 

s 

b 

7 

6 

u 

ill 

11 

12 

13 

14 

15 
lb 
17 
1 1 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 
5L 

51 

52 
c  3 
14 
5  5 

56 

57 
5P 
59 

0 
t  1 

f  Z 
t  3 


. COOOOQCOOl 1 1111  1111222222222233333  33333444444444455! 
01 2345676  901 2 345 6789C 123456 78901 2 34  567590123456789012 

DARK 


Ifi 


55556666 

67890123 

BRIGHT 

4 


*  *• 

*  **  *  * 

*  *  *****4*  *  *  *  * 

*  ******  *  **«****«•  *  *****  *  * 

***  *  ******************  *****  *  ****  *  *  *  * 
*  *  •  ***  **********  **********  *  **  *  * 

•  *•  •  ***************************  * 

A*****************************  *  *  *  * 

*  ****  «  ****** /\***********************  *  *  * 

*  *******  ***«**«A ****«& *******  A* A**0********  *•*  *  * 

*******A***A***AAAAA*A*«**AAAAAAAAAAA« ****** ********  * 

*  ********  AA**  A  AAA  AABBAPAAAABBC93CAA3AA***************  * 
»*  ******* a*BAB A AA8BA A AQC3CCB ACbCBBEBAA AA*A***********  *  * 
*******  A** AABBABBABCCCCCEDEODBOCOAABB A3  A ABAA*A A *4* ******** 

****A»A*A ABA A98BBB9SCC0EELEDDDC3CBCCBBBBACMBC** ******** A 
*****A*AA»A83CFDEDDrjDD£FEEDDDEDDDECDCDB3ACFQDB*AA**AAAAAE 

**  ***AABBCECCCCCCCEEFF0EDC0C0DC0CCDDCCCBDCCCB93BAAAABAAE 

*  ****** AABBBBCD0EED0BDC0DCCCDCBBCCBCCBCCDCB3C33BAAAAAA3 
**  **AABBBCBDOCDC0DCDL10BEDCDDCCCCDEDDCCEDEDDBCCCCC3AB*A J 

****AAABCCBC8DC0DBDEEEDEnDDDD:D0D0D3CEFEDDDDDD0EEE9CC3S 

•  **aaboocdccddefefedefeeeddedofeccededddefefdeodeeccecs 

-  *  ***ABCCEECDDEDEEELEO£DODDEOCCDDOOCCDCCEOCFCDCECEDEDDCCDDT 

**AABBOCCBCCCDBCAOOCDBCOCDDDBE3COnCCOOCBCCCB8BCCCC33B98W 

***aaacoecoddcodcbdeododoccccoddcdefdefefddeeededfecfeecds 

******* ABCnrDFDEOEDCFCCDCOCEFECDOEODDFEGFFEFGEFFFGGFFGBSOEO* 
***# AAACCCCCCCDDCCDCCCUCC00CD00CDC0CCDC0ECCCBCBCCC0DCDDC33S 
***********  AAA  A  APB  ABBB83  AAA  ACS  CB3t)ACBABABBBBB9BABB3BA**A  A  AL 

**  *******AAAAAAAAABAAAAAB6AA3BBAAA AABBA*AA ABA»*AAA*****A*AG 
.  *  ****AAAAA*B*AAAAB AAAAA AAAAAA AAAA AAA»A« AAAA**A**««A******F 
*****A**A*A***A**** A**A**AA*A**********A************A*****G 
******** *A*AAA***4* ***************************** **********3 
.  *  «**m**  ******  **********************************4*******C 
*  **  ************ *«A ******* *****«*************•**$ *******  4 

**  ********A****A **********  ********  *******  ********4 

******************** ****•«$***«  *********  ******  **  «******«3 

*  **  ******************  **  «***  ****  ***  **•*  **********  A 

**  **  ***  **  *****************  **  *******  *  ******  *  * 
!  **  ***  *****  **  ****  **  ******  *  ***  *****  *******  *  * 
*  ****  *******  ***  ***********  ****  ***  *** 


** 


*** 

*** 


**** 
•  *** 


** 


**  * 
** 


*  * 


WHERE  BLANKS:  0,  *=  I  To  9,  A=  10 
2-  260  TO  269,  AND  %-  270* 


TO  19,  3r  20  TO  29, 
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.19.  Cloud/clear/8ncw  decision  matrix. 
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Reliability  of  cloud/no  cloud  decisions  “  92. 6% 

()  *  Raw  frequency 

Table  3.1.  Verification  rates  using  the  decision  matrix  on  data  from  orbital 
revolution  2552  and  manually  prepared  analyses  as  "truth". 

MANUAL 

CLEAR  SNOW  CLOUD 

CLEAR  95.9%  0.2%  3.8X  (2294) 

S  SNOW  0.6%  94.9%  4.8%  (771) 

C  CLOUD  4.6%  4.6%  90.9%  (8212) 

(2580)  (1114)  (7586)  (11277) 

Reliability  of  cloud/no  cloud  decisions  92.3% 

()  ■  Raw  frequency 

Table  3.2.  Verification  rates  using  the  decision  matrix  on  data  from  orbital 
revolution  2680  (North  America)  and  manually  prepared  analyses  as  "truth". 
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()  *  Raw  frequency 


83.02 


Verification  rates  using  the  decision  matrix  on  data  from  orbital 
2860  and  manually  prepared  analyses  as  “truth". 
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Reliability  of  cloud/no  cloud  decisions  -  88.82 
()  *  Raw  frequency 

Verification  rates  using  the  decision  matrix  on  data  from  orbital 
2680  (Antarctica)  and  manually  prepared  analyocs  .•’8  "truth". 
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Table  3.5.  Verification  rates  using  the  decision  matrix  on  data  frca  all  or¬ 
bital  revolutions  and  manually  prepared  analyses  as  "truth". 
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contain*  normalised  distribution*  of  the  result*  froa  the  algoritha  versus  e 
given  verified  result.  The  rata  of  success  (reliability)  of  cloud/no  cloud 
decisions  is  shown  at  the  bottom  of  each  figure.  When  we  computed  these 
reliabilities,  we  considered  clear  and  snow  categories  sa  cloud-free  because 
neither  one  contained  cloud  within  tha  IFGV.  The  only  difference  was  that  one 
contained  snow  and  tha  other  one  dii  not.  The  overall  reliability  was  an 
excellent  90.12. 
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SECTION  4  -  LIMITATIONS 

While  completing  this  study,  we  noticed  several  factors  which  limited  the 
utility  of  SSC  data.  Limitations  to  using  SSC  data  include  sensor 
characteristics,  data  flow  difficulties,  and  inherent  technique  limitations. 
The  first  two  classes  of  limitations  can  be  removed  or  reduced  by  appropriate 
engineering  modifications  to  follow-on  sensors.  Below  ve  discuss  three 
aspects  of  the  limitations  due  to  sensor  characteristics.  Then,  we  address 
the  data  flow  difficulties  and  irherent  technique  limitations. 

We  estimated  that  approximately  10Z  of  the  SSC  samples  which  were 
categorized  as  cloudy  were  only  partially  filled  with  cloud.  In  some  cases, 
only  one  L-detector  pixel  within  the  IFOV  contained  cloud.  We  learned  from 
experience  that  L-detector  data  were  much  sharper  than  SSC  data.  If  the  SSC 
and  L-detectors  had  the  same  resolution,  there  would  have  been  less  overlap  of 
SSC  cloud  samples  into  the  clear  and  snow  categories.  We  believe  this  would 
improve  the  reliability  to  95%. 

During  pre-launch  preparations,  WEC  calibrated  only  two  of  the  48  detector 
elements  of  the  SSC.  The  resultant  inter-channel  response  to  a  given  scene 
radiance  varied  by  an  estimated  +5Z.  This  inter-channel  variability  was  small 
enough  when  compared  to  the  wide  separation  between  joint  L-detector  and  SSC 
responses  from  cloud,  snow  (ice),  and  clear  scenes  that  the  affect  on 
reliability  was  negligible. 

The  design  of  the  SSC  gain  control  system  did  not  provide  the  capability 
to  reduce  the  sensitivity  of  the  sensor  to  eliminate  the  saturation  ceu»^4 
cloud-free,  arid  terrain  features.  Accordingly,  we  could  not  prove  ov.r 
hypothesis  that  the  linear  relationship  between  SSC  anJ  L-detector  responses 
also  applies  to  theee  features.  As  a  result,  our  cloud/no  cloud  decision 
algorithm  erroneously  labeled  more  clear,  arid  features  aa  cloudy  at  SSC 
saturation,  and  vice  versa,  than  would  have  been  the  case  with  better  gain 
control. 

Occasionally  some  SSC  samples  were  lost  during  playback  and  transmission 
due  to  noise  in  the  recorded  data.  Unfortunately,  we  could  not  identify  which 
scan  lines  had  been  lost  because  data  time  information  was  only  provided  once 
each  minute.  As  a  result  some  SSC  and  L-detector  data  were  mismatched. 
However,  we  were  careful  to  select  cases  for  verification  where  this  effect 
was  minimized. 

The  SSC  was  designed  to  measure  reflected  solar,  near-IR  energy.  As  a 
result,  we  used  data  from  daylight  areas  only  in  this  study.  We  were  not  able 
to  obtain  meaningful  results  in  areas  of  twilight  or  total  darkness  (i.e., 
areas  with  zero,  or  less,  solar  elevation  angle.) 
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SECTION  5  -  CONCLUSIONS 


The  germanium  detectors  were  stable  throughout  nearly  7  months  of  on-orbit 
operation.  No  resolvable  deterioration  in  the  performance  of  the  sensor 
occurred  during  this  period. 

Signal  to  noise  ratios  were  high  enough  to  permit  processing  of  matching 
SSC  and  L-datector  samples  in  daylight  from  near  solar  zenith  to  the 
terminator. 

Data  from  a  narrow-band  infrared-channel  centered  near  1.5um,  used  in 
conjunction  with  L-detector  data,  will  detect  clouds  with  more  than  90X 
reliability. 
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SECTION  6  -  POTENTIAL  BENEFITS 


The  second  objective  of  our  study  was  to  identify  potential  benefits  of 
processing  snow/cloud  discriminator  data.  Specifically,  we  looked  for 
techniques  and  data  sources  to  improve  the  Three  Dimensional  Nephana lysis 
Model  (3DNEPH)  and  the  Snow  Cover  Analysis  Model  (SNODEP).  We  have  listed  the 
most  likely  benefits  in  the  following  paragraphs. 

The  greatest  improvement  in  3DNEPH  analyses  would  result  from  better 
detection  of  low  stratus  clouds  over  snow.  Usually,  the  3DNEPH  algorithm 
which  analyses  infrared  data  does  not  detect  these  clouds  because  the 
temperature  of  the  cloud  tops  is  frequently  equal  to  or  warmer  than  the 
temperature  of  the  underlying  snow.  Figure  3.12  contains  a  picture  of  nearly 
26,000nm2  of  stratus  clouds  which  the  3DNEPH  algorithm  would  have  been 
unable  to  detect.  However,  our  cloud/no  cloud  decision  would  have  detected 
these  clouds  with  a  reliability  of  better  than  902. 

The  greatest  improvement  in  SNODEP  analyses  would  result  from  the  use  of 
additional  snow  cover  information  in  areas  of  the  world  where  there  are  no 
stations  reporting  snow  depth.  The  snow/no  snow  decision  algorithm  could 
provide  snow  cover  information  within  these  areas  in  lieu  of  the 
climatological  data  now  used  by  SNODEP. 
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AWS  DISTRIBUTION 


DISTRIBUTION: 

1WW  (10) 

2WW  (10) 

3WW  (  6) 

5WW  (  7) 

7WW  (  5) 

AWS/DNT  (1) 

2WS  (1) 

AFGWC  (5) 
USAFETAC/DN  (2) 
USAFETAC/OIi-A  (1) 
USAFETAC/TSK  (5) 

AUL  (1) 

3350  TCHTG/TTMV  (5) 
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